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Abstract 

We define the notion of A-model Lagrangian D-branes as introducing defects in the 
Calabi-Yau crystal. The crystal melting in the presence of these defects reproduces all 
genus string amplitudes as well as leads to additional non-perturbative terms. 
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1 Introduction 

In [1] a duality between topological A-type string theory and statistical model of a three 
dimensional melting crystal was proposed. In particular it was demonstrated that partition 
function Z crysta i of the melting crystal corner computes the closed string partition function 
Z c i ose d in the C 3 background. More precisely, 

Zcrystai(q) = Z dosed (q), q = exp(-g s ) (1) 

where the string coupling constant g s is identified with the inverse of the temperature of the 

crystal. Also in [1] it was explained how to compute the topological vertex [2] in the crystal 

picture. In [4] this duality was interpreted as coming from quantum gravitational foam 

for target description of the A-model strings, which is a "quantum Kahler gravity". Each 

geometry contributing to the quantum foam corresponds to a configuration of molten crystal 

(roughly speaking each atom removal from the crystal changes the quantum foam by blowing 

up of a corresponding point). It was proposed in [4] that in the limit j s <C 1 the quantum 

i 

foam (the melting crystal ) is a fundamental description of physics at short distances ~ gih- 



1 



String emerges in the high temperature limit (1/T = g s « 1) as an effective description of 
the theory at the string scale l s . A related mathematical proposal was also presented in [5]. 
Further developments have shed light on this duality: It was conjectured in [6] that A and 
B models on the same Calabi-Yau are S-dual to one another. Moreover it was shown in [7] 
that this S-duality follows from type IIB S-duality and that the atoms of the crystal can also 
be viewed as D(-l) instantons of the B-model. See also [9]. 

In this note we continue to investigate the proposed duality and would like to give the 
statistical interpretation for non-compact Lagrangian D-branes in the context of C 3 . It turns 
out that the statistical interpretation of Lagrangian D-branes is introducing defects to the 
crystal. For the case of branes at some particular framing we find that the topological string 
amplitudes are identical to that of melting crystal in the presence of the defects introduced 
by D-branes. We also find the defect interpretation of the lattice for arbitrary framing and 
also in the presence of anti-D-branes. However for these cases the corresponding melting 
rule needs to be modified, similar to the melting rule in the context of local toric 3-folds 
(when we have more than one corner for the lattice). We will make some comments about 
what this melting rule would look like for these more general branes, but will not develop it 
in full detail in this paper. 

The organization of this paper is as follows: In section 2 we review the results of [1]. In 
section 3 we present the geometry of modified crystal for branes and anti-branes and also 
give melting crystal picture for certain configurations of D-branes and anti-D-branes. In 
section 4 we discuss general brane and anti-brane configurations at arbitrary framing. In 
section 5 we show that our melting crystal picture of the D-branes is compatible with the 
semiclassical limit of small string coupling. 

2 Review of the closed topological strings /melting 
crystal duality 

In this section we review the duality between closed A-model topological string theory on 
C 3 and statistical model of melting crystal corner proposed in [1]. The A-model topological 
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string theory studies holomorphic maps from Riemann surfaces to the Kahler target space. 

$ : E -> X (2) 

The partition function of this theory is given by 

Z A = exp{Y.9 2 r 2 ^), T g =\ e~ A ° (3) 

where A c is the area of the holomorphic curve C = 3>(£ g ) of genus g, and integration goes 
over the moduli space of maps Ai map .s,g with a suitable measure. Here g s denotes the string 
coupling constant. When the volume of X is large the leading contribution to Za comes 
from constant maps and the moduli space becomes a product M. ma p S ,g = M. g x X, where 
M. g is the moduli space of genus g curves. In this large volume limit the expressions for T g 
were obtained in [11,12]. As follows from the results of [11, 12] the closed string partition 
function on C 3 (taking its Euler character to be 2) is given by the McMahon function: 

oo 

Z c i osed {q) = M(q) := II (! " O"", Q = ^ ( 4 ) 

n=l 

Now let us recall how the same function arises in the dual melting crystal picture [1]. Consider 
cubic lattice in the positive octant of R s with lattice spacing given by g s and define the rules 
of melting. One can remove atom with coordinates (xo, yo, zq) from the lattice if and only if 
all the following sites are already vacant: 

(x,y ,z ) x<x , (x ,y,z ) y < y , (x ,y ,z) z < z (5) 

With these rules, each configuration of k molten atoms corresponds to a "3d Young diagram" 
(a 3d partition) of k boxes in the positive octant of R 3 . Now we sum over the random 3d 
partitions with the weight q# boxes ; where q = e~ 9s and g s stands for inverse temperature of 
the crystal (measured in units of chemical potential for removing atoms). To compute 

Z crys tal{q)= E ( 6 ) 

3d partition n 

we recall that, as discussed in [10], a 3-dimensional partition n can be thought as a sequence 
of 2d partitions {fi(t)}, t e Z obtained from diagonal slicing of n by planes x — y — t. These 
2d partitions obey the interlacing condition: 

fi(t) < fi(t + 1), t<0, fi(t + 1) < fi(t) t>0 (7) 
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where the 2d partitions \x and v interlace, and we write [i > u, if 

A*i > f 1 > A*2 > v 2 ■ ■ ■ (8) 

Here and below we let fik denote the number of squares in the k-th row of //. The diagonal 
slicing allows us to use a well known map from 2d partitions to states in the NS sector of a 
complex fermion: 



^w = Et + ^" n " 1 . ff^Ec+i^ 1 . {Vw±>v>: (m+ i)} = <w (9) 

nez neZ 

The state corresponding to a partition \i is given by 

i^n^iMjo) (io) 



1=1 

where d is the number of squares on the diagonal of [i and, denoting the transposed 2d 
partition by ji T , we define 

di = Hi-i + ^, bi = /iJ-i + ~ (11) 

The example of the state/2d partition correspondence is shown in Figure H 
Now, Z crysta i(q) can be easily recast in the operator language: 



Zcrystalyl) — ( 



n<^r + (i)W n r_(i) g 



L 

1 -A-s/y 

i=0 ' x i=-oo 



(12) 



Here |0 > is the state annihilated by all positive modes of ip and ip* and operators L , r±(l) 
act in the following way. 

LolfA >= \fi > (13) 



r4i)\p>=52\u>, r + (i)\n>=J2W> (u) 

where stands for the number of squares in 2d partition \i. The final step in computing 
Z cr ystai(<l) is to rewrite the correlator (|T2*|) in the form 

oo oo 

z cryst ai(q) =< o| n r+ n r_(g- m+ 5)|o > (15) 

n=l m=l 




Figure 1: We draw 2d partition fi in a diagonal slice at t — —2. The numbers of rows are 
Hi =4, /i 2 = 2, /i 3 = 1. The fermion state is |/i >= ^-7/2^— 5/2 / 0^i/2'V ; — 1/2 10 > . 



where 



T + (<f-*) = q-^ Lo T + (l)q {n -^ u \ V -{q' n+ ^) = q-^+^T 4l)q^ n+ ^ )L ° (16) 

We will say more about the operators T±(z) and their simple realization in the bosonized 
picture in the next section. Here we only need their commutation relation 



r + (z)r^z') = (i-z/z'y 1 T4z')r + (z) 

which leads immediately to : 

Z cry staM) = X{(l-q n Y\ q = e~ B ' 

n=l 



(17) 



;is) 



Comparing (|3J) with (|18jl finishes the proof of the duality (pQ) for closed strings. In [1] it 
was shown how one can formulate the topological vertex amplitudes [2] in terms of crystal 
melting using the above operator formalism. 
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Figure 2: Projection to the base of a D-brane ending on y-axis at y=a. 

3 String/Crystal duality with D-branes 

In this section we give the melting crystal interpretation of D-branes in topological A- 
type string theory on C 3 . We first recall the relevant results from string theory. In particular 
we explain how the D-branes affect the integral of Kahler form around 2-cycles surrounding 
D-branes and we use this to formulate the D-branes in the crystal melting problem as in- 
troducing defects. We show that the amplitudes for D-branes agree to all orders in string 
perturbation with the melting crystal picture. We also get additional non-perturbative 
pieces, suggesting a non-perturbative completion of string theory amplitudes. 

3.1 D-branes in topological A-model string theory on C 3 

One may view C 3 as a T 3 fibration. Let (zi,z 2 ,z 3 ) denote complex coordinates of C 3 
and let 

(x,y,z) = Oi| 2 , \z 2 \ 2 , \z 3 \ 2 ) 

be the coordinates on the base + which is the positive octant in R 3 . We denote the coor- 
dinates on the fiber T 3 by 8i, i = 1, 2, 3. Now consider a Lagrangian D-brane which has the 
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geometry given by 

y = x + a = z + a, 6»i + 2 + 3 = 0, a > 0. (19) 

The D-brane has topology S* 1 x i? 2 and ends on the y-axis at y — a ( see Figure EJ). As 
discussed in [2] there is back-reaction of the geometry to the presence of the D-brane. Namely, 
the periods of the Kahler form K are changed by the amount 

Aj^K = g s (20) 

where £ is any 2-cycle linking the D-brane world-volume. This back-reaction is going to be 
our guide in search for crystal interpretation of the brane. 

Also, as explained in [14], in order to define quantum theory on a non-compact D-brane 
one has to specify boundary conditions at infinity called the "framing" of the D-brane. The 
trick is to modify the background by introducing extra direction / in the toric base over 
which one of the cycles of the T 3 fiber degenerates along a "brane" in that direction at 
infinity. There is an integral ambiguity in the choice of /: 

p is called the framing number. This leads to the notion of a preferred plane containing 
the D-brane. Let's call this plane the D-brane plane Pp. A choice of the framing is also 
equivalent to the choice of the plane Pd containing the D-brane and the choice of one of the 
two directions mi,rri2- This follows from the orientation of / versus — /. The concept of 
framing is illustrated in Figure HU We will show that the notion of framing enters naturally 
in the crystal melting picture of the D-brane. 

The choice of p = framing in the convention of [14] 1 is the simplest one. The partition 
function of the topological A-model string theory on C 3 in the presence of the D-brane (anti 
D-brane) with geometry (|T9|) and framing p = is given by [13] 

OO , N 

Z% Hn9 (q, a) = M(q) H(l- e~ a q n ^ ) , q = e"* (21) 
1 p = —1 in the convention of [2] 
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Z§ rm9 (q, a) = M(q) e~ a q^ ) (22) 



n=l 



where M(q) is the McMahon function (£Q). Note that since the brane amplitude at framing 
p is equivalent to anti-brane amplitude at framing 1 — p the above equations can also be 
viewed as statement for anti-brane (brane) with framing number 1. 

Also, note that free energy F = —logZ^ rm9 (q, a) has the following form: 

co —na 

F=T T7 TT (23) 

^! 2nsinh(n# s /2) v ; 

so that perturbative expansion F g ~ (g s /a) 29 ~ l breaks down for a ~ g s . Nevertheless, the 
partition function Z^ rm9 (q } a) is well defined for a ~ g s and it is natural to ask what this 
sum computes. To find the answer to this question we are going to study the interpretation 
of D-branes in the melting crystal. 

3.2 D-branes as defects in the melting crystal 

How can one detect the presence of a D-brane in the crystal melting picture? To answer 
this question we would like to use the important information about the back-reaction of the 
geometry (|2(J|). Let us first recall the role of the Kahler form K in the crystal picture. It was 
discovered in [4] that the melting crystal partition function Z cryst ai{q) arises from the sum 
over generalized 2 Kahler geometries on C 3 with Kahler forms K quantized in units of g s . 
From this viewpoint the number of atoms removed from the corner can be thought of as the 
change in the number of holomorphic sections of the U(l) bundle whose curvature is K/g s . 

Now we introduce a D-brane of geometry (jl9j) . The projection of the brane to the base 
and the framing vector define a plane Pp in + . This plane intersects z = plane along 
the half-line mi : y = a + (1 + p)x and x = plane along the half-line m2 : y = a — pz 
(see Figure EI). There is a natural choice of a 2-cycle linking this D-brane which we denote 
by E dictated by the choice of framing: This cycle is a difference of planes whose projection 
to the toric base gives a pair of half-lines Lj, L 2 parallel to Pjj in the direction of mi (fixed 
by the choice of framing). We postulate that this gives rise to a 'closed cycle' at infinity 

2 In a sense that the sum includes non-geometric excitations required for a consistent quantum theory 



z 



X 




Figure 4: The half-lines Li and L 2 enter the definition of the 2-cycle E. A cross stands for 
vacancy and dot for atom. 

E = L 2 — Li, compatible with the fact that at infinity the framing of the D-brane chooses 
an extra vanishing direction of a circle. Here L\, L 2 are given by 



where 7/1,7/2 and z* are fixed nonnegative integers such that y\ < a < y 2 . The lines L\ and 
L 2 together with 81 angle give rise to planes, which for general p are not holomorphically 
embedded in C 3 . The period of the Kaler form K along the 2-cycle E is given by 



where \Li\ is the number of atoms on the half-line Li, i — 1, 2. From (J23|) follows that the 
change of the period by one in g s units, which is the characteristic signature of the D-brane, 
can be achieved by, for example, dropping a single atom from the half-line L\. We illustrate 
this idea in Figure 0] 

Recalling that parameters yi,y 2 ,z* which specify a linking 2-cycle E are arbitrary non- 
negative integers ( with the only constraint y\ < a < y 2 ) we see that in the presence of a 
D-brane one layer of atoms is deleted from the x = plane, in the region < y < a. In 



L\ = (x,y 1 + (l+p)x,z*) 



L 2 = (x,y 2 + (l+p)x,z*) 



(24) 




(25) 
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Figure 5: The shape of the melting crystal in the presence of D-branes. The arrows carry 
negative (positive) sign and show the evolution with T_ (T + ). The crosses are vacant sites 
inside the cylinder Ai u . 

other words the crystal develops an extra corner. Thus, we are led to the interpretation of 
D-branes as defects modifying the shape of the crystal. For many D-branes the crystal will 
develop many corners. Moreover the concavity or convexity of the corner is correlated with 
whether we have a D-brane or anti-D-brane. See Figure 

Note that the picture of modified crystal with D-branes at various framings makes sense 
as long as the corresponding m 2 lines do not intersect the crystal axis z or the other corners 
developed. Thus only for these cases we can propose a simple dual crystal interpretation. 
Moreover the notion of linking cycle £ that we defined will only make sense for two choices 
of framings for each brane (p = 0, —1) otherwise the corresponding cycles intersect the 
boundary of the crystal. This already suggests that perhaps the rule of crystal melting may 
be different for general framings and this indeed is consistent with what we shall find. In 
fact the simplest rule is for framing p = and that is the case we focus on. 

One would naturally propose that the modification due to the D-brane in the A-model 
partition function is the same as the statitical mechanical model of melting crystal but now 
with the modified crystal + /A4 u studied in [1,10]. Here Ai u is a cylinder with the 2d 
partition v as its base. In other words, one considers 3d partitions in + and regards the 
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two 3d partitions as identical if they differ only inside M. v . The weight of the 3d partition is 
determined by the number of boxes outside M. v . We will see that this is indeed the correct 
rule for the case of framing p = 0. For general framing the crystal geometry is modified as 
we have found, but the melting rules also need to be modified. Some aspects of this will be 
discussed in the next section. In the remainder of this section we concentrate on the case 
with p = and show how to recover the D-brane amplitudes from the melting crystal. 

In the operator language the 3d partition in the modified container is obtained from 
alternating evolution with T_ and T + whose action on diagonal slices provide appropriate 
interlacing pattern. We illustrate the evolution giving "room with the corners " in Figure 
[S| Now we would like to check the proposed interpretation for the D-branes by explicit 
computation using realization of F± as operators in fermion Hilbert space. It turns out 
that this realization of D-branes is consistent with the intuition developed in the B-model 
topological string theory [2] that non-compact D-branes are fermions. 

Let us recall [10] the representation of operators T± in the bosonized picture. One 
introduces chiral boson <p(z) such that 



P(z) 



=: e 



z 



(26) 



and expand 4>(z) into zero, positive and negative modes 



<t>(z) = 0oO) +0+0) +<f>-(z), <f>+(z) } <f>-(w) =log(l 



z 



(27) 



w 



Then T±(z) are given by 




(28) 



Note that ip(z) and ip*{z) are related to T±(z) : 



= e-^r_(^)r; 1 (z), i/j(z) = e^VZ 



(29) 



Now we propose that introducing the D-brane of the geometry (fT9"j) at a = g s (N + |) with 
framing p=0 amounts to inserting the operator 



* D (z)=rz 1 (z)r+( z ) 



(30) 
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at z = q~( No+ h) in the correlator (|T3|) . We choose not to include the zero mode part 
in the definition of the D-brane operator ^ d{ z ) since from the crystal viewpoint it is more 
natural to deal with zero fermion number correlators. As discussed in [8] disregarding zero 
modes corresponds to having compensating (fermion number ) flux at infinity. 
To verify our proposal for the D-brane operator let us compute 

oo No + 1 oo 

z^ ta/ ( g ,iv ):=<oinr + (g n ^) n r_(r m+ ^(g- (iVo+ ^) [] r_( g — +l)|o> 

n=l m=l ' N +2 

(31) 

Using our definition (|50|l we have 

oo No oo 

z c ^ stal ( q ,N )=<o\l[r + ( q ^) nr_(?- ro ^)r + (fM)) j| r_( ff -™ + i)|o> (32) 

n=l m=l ' ^0+2 

Note that insertion of a D-brane creates an anticipated disorder at q~~( No+ k) in a sequence 
of T_ operators. The commutation relation (fTTj) gives the result 

OO / s 

ZZ ysta \q,N G ) = M(q)aq) U U - e^ ^)^ , g = e"* (33) 

n=l ^ ' 

where £(g) = n^=i(l ~~ g™) -1 - Let us compare string answer (j2~Tj) with (J3~3~j) : 

^ ta/ (g, No) = mZ% rm9 (q, a), a = g s (N + i) (34) 

The renormalization factor £ (g) is important for the melting crystal picture of the D-brane. 
Due to £(g) factor we can write Z^ ystal (q, No) in the simple form: 

No 

Z^ stal (q } No)=M(q) J] (1 - O" 1 = P(o, 0, i/(JV )) (?) (35) 

n=l 

Here v(No) stands for the 2d partition with single row of length JV and P(o, 0, uj (q) was 
defined in [1] as the partition function of the 3d partitions in the modified container + / M. v . 

P(o, 0, u(Nofj (q) can also be viewed as the partition function of the melting crystal when 
melting starts from the cubic lattice A G + with all the sites 

(Q,y,z), y < N o9s 

already vacant (see Figure |H1 ) • 
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Figure 6: Modified lattice for a D-brane with p = 0, N = 3. The crosses stand for vacancies 
and dots for atoms. 



We now easily generalize the above result to give melting crystal picture for the M D- 
branes of the geometry (fElj) with a* = g s (Ni + |), i = 1, . . . , M and Nm < Nm-i < ■ ■ ■ N±. 
In the operator language we compute 

oo N M +l 

ZZTLnes^^) :=< 0| J[ T + (q^) R T D {q^^) (36) 

n=l m=l 

IJ r_(g- m+ ^ D (V^M-i+^) rj r_( g - m+ *)... 

m=N M +2 m=N M _ 1 +2 
Ni+1 oo 

... f[ r_(g- m+ 5)^(g-W+5)) f[ r_(g- m+ 5)|0> 

m=N 2 +2 m=N 1 +2 

Evaluating the above correlator using (JT7I ) we find 

z^n-LnesiQ, Ni) = M( q ) n i - g^-^ n n (i - e"*)- 1 = Wo, o, ^ (3?) 

i<j ^ ' i=l ni=l V ' 

where v^N^ stands for the 2d partition such that z/j = iV, — M + z. The last equality in ()37p 
follows from the relation established in [1] 

P(0, 0, !/) (q) = M(q)q-^C. M (q- 1 ) (38) 
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Figure 7: Modified lattice for two D-brane with p = and N\ = 6, N2 = 3. 

where ||A|| 2 = J2i an d C..!/(<7 -1 ) is the topological vertex defined in [3]. In the crystal 
language one is melting the cubic lattice in the positive octant with all the sites 

(g a (i-l),y,z), y<g s (N i -M + i), i = l,...,M 

already vacant (see Figure [7| for the case of two D branes) . 

We have considered D-branes which end on the y-axis. Let us now insert an anti-D-brane 
ending at x — g s (No + |). We claim that the anti-D-brane operator is given by 

= r-^r; 1 ^) (39) 

inserted at z = q( No+ ^ where subscript x indicates that we are dealing with anti-D-branes 
ending at x-axis. This is clear from the jumping rule for the Kahler class. We can now check 
this claim by computing the amplitude: 

00 -/Vq+1 00 

Z^ sW (q,N Q ) :=< 0| II r + (g^)^ (<?<*>+*>) r+(<T^) U r_(<T n+ *)|0 > 

m=No+2 m=l n=l 

(40) 

We use (f3~H|) to recast Z^ y * tal (q, N ) as 

00 No 00 

z^ stal (q,N ) =< o| n r + ((T-i)r_( g w+i)) n r+(<T^) n r_(g-" + *)|o > (4i) 

m=Na+2 m=l n=l 
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Figure 8: Modified lattice for anti-D-brane with p — 1, Nq — 2 ending at x-axis. 

Note that insertion of x creates a disorder at q( No+ ^ in a sequence of T + operators. 
Computing (j41)l we find 

Zgf a '(g,iVo) = II 1 " e~"'^ + ^A q = e"* (42) 

n=l ^ ' 

where = n^=i(l ~~ g n ) _1 . Note that is equal ( up to the factor £(g) ) to the string 
answer (|21j) for anti-D-brane with framing p = 1. 

Z^ y ^ tal (q, Nq) can also be viewed as the partition function of the melting crystal when 
melting starts from the lattice A x £ (9 + with all the sites 

(x, 0,2;), x < N g s 

already vacant (see Figure |HJ). 

We would like to end this section with a comment about extra factor £ (g) which appears 
in the crystal melting picture as compared with the string theory answer (|34|). To ensure 
that this extra factor cannot be detected in perturbative stringy computations we propose 
to modify the definition of a D-brane operator in crystal picture 

*d = q- 1/2 ^D (43) 
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Then, extra factor would be 77(g) 1 and the change in the free energy due to this factor is 
non-perturbative in the g s — > limit. 

AF = -logrj(q) = -logrj(q) + ~log^-, q = exp(-^L) (44) 

Z Z7T v g s ' 

In the S-dual language [7] the factor 77(g) -1 may be thought as coming from extra D(— 1) 
instantons in the presence of Lagrangian NS2-brane. It would be interesting to determine 
analogs of £(q) for arbitrary framings. 



4 Crystal amplitudes for more general D-brane config- 
urations 

So far we have discussed crystal melting interpretation for special configuration of D- 
branes. Here we would like to consider more general configuration of branes and anti-branes 
at various framings. 

4.1 Anti-D-brane amplitudes 

Let us consider anti-D-brane ending on the y axis and zero framing. We will take into 
consideration that the jump of the Kahler period along the 2-cycle linking anti-D-brane is 
opposite in sign to the corresponding jump around D-brane. Then, we can naturally propose 
the shape of the melting crystal as shown in Figure 121 Note that the corners are now concave, 
as opposed to the case for branes ending on the y-axis. 

To verify this proposal and find the rules of melting we again use operator language. We 
suggest that anti-D-brane ending on y-axis at y — g s (No + |) is described by insertion of the 
operator 

^(z)=T_(z)T- 1 (z) (45) 

at z = q~( N °+h\ 

OO ;Y() — 1 OO 

Z ^\q } N ):=<Q\Y[T + {q^) J] T_(q- m ^)^(q-^) r_(g- + l)|0> 

n=l m=l N +2 

(46) 
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z 








Figure 9: The shape of the melting crystal in presence of anti-D-branes. 
We use (|43|) to recast Z^ ystal (q, N ) as 

oo No oo 

< o| n r +(9 n_l ) II r_(g- m+ 5)r 2 _(g-W' + 5))r; 1 (g-( 7Vo+ 5)) f[ r_( g - m+ *)|o > (47) 

n=l m=l ' N +2 

Note that insertion of an anti-D-brane tyjj creates a much more complicated disorder at 
q~( No+ k^ in a sequence of T_ operators than D-brane does. The commutation relation (|17j) 
gives the result for anti-D-brane 

oo , \ -1 N 

Z^ sta \q,N Q ) = M{q)C\q)\{[l-e-^ N ^h n -A = M{q) ]J (l-g»), g = e ~* (48) 

n=l ^ ' n=l 

Let us compare string answer ({22} with (j48jl : 

at ) = r i( g )z£**( g> a), a = g s (N + i) (49) 

Now we would like to give a crystal interpretation of (|48j). A natural guess which gives the 
correct jump of the Kahler class is to start melting the cubic lattice A in the positive octant 
with vacant sites 

(0,y,z), y>g s N . 

It turns out that in addition to the choice of this modified lattice, one also has to modify the 
melting rule to reproduce the correct answer predicted from the operator formulation (and 
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known string amplitudes). To formulate the modified rule we use 



II(i-g n )= E (-)^=>^=i m * 



(50) 



n=l ^e"Pjv, 



where Ptv is a set of column strict 3 2d partitions with maximum column length l max < Nq 
and rii = 0, 1 counts if a column of length i is present in Here is a tentative proposal for the 
modified rule. The melting in x = plane starts from ^0, (No — l)g s , and is going in the 
direction of decreasing y and increasing z. The allowed configurations for this new melting 
are collections of atoms forming column strict 2d partitions. So for example, a configuration 



is allowed. The sign of configuration is the number of columns, which is two in the example 

()51|). We give a more nontrivial example of an allowed configuration molten in the x = 

plane in Figure The melting in the rest of the positive octant goes in a standard way 

© independently of the melting in the x = plane. For this reason a total configuration 

of molten atoms has the form (A; B) where A is a configuration obtained from melting in 

+ and B is a configuration (with signs) in the x = plane. For more anti-D-branes or 

mixtures of D-branes and anti-D-branes, one naturally expects that the modified lattice be 

consistent with the change in Kahler class. Viewing D-brane as fermion and anti-D-brane 

as anti-fermion, this is also consistent with the fact that j> d(f) measures the Kahler class of 

the crystal and its jump across the insertion of fermion or anti-fermion is consistent with the 

expected geometry of the lattice. However, the precise rule for the melting crystal remains 

to be worked out for this more general melting. Of course the operator expressions always 

exist for such amplitudes. In this regard, we would like to emphasize that the global rule 

of melting is also absent for closed strings on toric Calabi-Yau manifolds with more than 

one fixed point of the toric action. The rules given in [4] are local 4 near each of the fixed 
3 Column strict means that the lengths of columns are strictly decreasing 

4 There are also modifications which account for world sheet instantons wrapped on compact 2-cycles 
between fixed points 




is not allowed but a configuration 




(51) 
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Figure 10: An example of an allowed configuration molten in x = plane. The number of 
columns is three, N = 8, the contribution to Z is (— l) 3 ^ 11 . 
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Figure 11: N% = 3, N2 = 11. 6 atoms are molten near each of the branes on top of the 
three molten strings. 

point. This is related to the fact that what appears as brane to one corner of the crystal will 
appear as anti-brane to another corner. In the spirit of [4] we can get a rough idea about 
the melting crystal interpretation of the D-brane and anti-D-brane when they are far apart. 
In this case we can regard our models of melting as local ones which are valid in the vicinity 
of each brane. 

Let us consider a D-brane at a = g s (N 1 + |) and an anti-D-brane at b = g s (N2 + 5) so 
that iV 2 > iVi. Note that in the presence of D-brane and anti-D brane there are essentially 
new types of objects (one dimensional lines) to be added to the melting picture [4]. These 
are "strings of atoms " of length 5 N 2 — N\ . The new configurations account for the factor 
5 the number of atoms in a string 
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(1 — q N ' 2 Nl ) 1 present in the partition function: 

N\ N 2 

Z D n = M(g)(l - J] (1 - g")- 1 J] (1 - g ™) (52) 

n=l n=l 

We show in the Figure ^2 a typical configuration of molten atoms on top of melting three 
string-like objects. 

4.2 The D-brane framing in crystal picture 

Here we would like to discuss the crystal amplitudes for D-branes in arbitrary framing. 
We will show, very much like the case for the anti-D-branes, that the operator formulation 
of the amplitude is compatible with the modified lattice we had obtained before. Though a 
direct interpretation of the statistical mechanical model encoded by the operator formulation 
is unclear. 

Let us recall the string theory partition function in the presence of a D-brane (anti-D- 
brane) of the geometry (j!9|) with framing p [2] 6 . 

oo k 

Z% rin9 (q,a,p) = M(q)Y,e' ka (-l) {1 - p)k q {1 - p)k{k - 1)/2 q k/2 Yl( 1 ~ Q = e" Ss (53) 

k=0 ra=l 

oo k 

Zi rm9 (q,a,p) = M(q)Y,e- ka (-lY k q^ k - l)/2 q k/2 \{^ ~ <?T\ (54) 

k=0 n=l 

Now we are going to illustrate how string answers can be reproduced from the operator 
formulation on the modified crystal (up to the renormalization factors £ ±x ). We define the 
fermion operators ip^ and ip*^ : 

^p)( z ) = i} n+ i(-) p{n+1) q~ Mn+1)/2 z~ n ~\ i> <P \z) = - ^ +i (-)P( n+1 )gP n ( n + 1 )/ 2 2-"- 1 

ngZ " nGZ 2 

(55) 

and in analogy with (j3U46|) 7 insert them at q~^ No+ ^ 

oo iVo+1 oo 

Z c ^ stal (q,N ,p) :=< Vl |nr + ( 9 ^) II r_(g- m+ ^)^)(g-(^)) [] r_(g— +i)|0 > 

n=l m=l ' ' N +2 

(56) 



6 We use convention for framing as in [14] 

7 It is more convenient here to insert operators with nonzero fermion number and thus change the vacuum 
at t = oo. For a single D-brane (anti-D-brane) this is equivalent to our treatment in section 3. 
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^■-D(^Ar ,p):=<v_ 1 inr + (g"^) n r.( q - m ^)r {p) (q' {No+ ^) n r_( g - + i)|o> 

n=l m=l N +2 

(57) 

where < 0| =< v^e* and < 0| =< v_ 1 |e _< ^°. In order to compute (JSTIj) we first find 
correlators involving modes ip n+ i'- 

oo No + 1 oo 

V n =< v x | n r + (g n ^) I] r_(<T m+l Vn + § II r_(<T m ^)|0 > (58) 

n=l m=l iV +2 

These can be extracted from the known expression for Z£ ystal (q, N 0l p = 0) © to be 

k 

V n= „ 1+k = M(q)aq)(-) k q k2/2 f[(l-q s r\ k > 0, P„ = 0, n < -2 (59) 

For anti-D-brane the procedure is analogous. Then we immediately find the relation between 
crystal and string partition functions 

Z- ysW (q, N ,p) = mZn Hn9 (q, a,p), a = g s (N + J) (60) 



Z^ stal (q,N ,p) = r 1 (q)Z^ rm9 (q,a,p), a = g s (N + I) (61) 

where £(q) = Yl™ = i(l — q n )~ l ■ Now an important point is that we can write operators ip^ 
and V* (p) as [2, 8] 

^ = u p *pu- p , r {p) = u P ru- P (62) 

Here U p = (-l)P L o gP c 2 and 

L = £ (fc + 1) : Vw^-*-A C 2 = - £ : ^ + lC fe _i : (63) 

fcez fcez z 

Recall that L and C*2 act on the state \fi > corresponding to 2d partition /i as follows: 



L \fJL >= \fi\\fJL >, C 2 \n >= -Kplfi > (64) 



where |/i| = J2i Hi an d = J2(i,j)e^^~J- As was discussed in [1] operator U_ p is responsible 
for the change of slicing of the 3d partition. More precisely, £/_ p maps a 2d partition in a 
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Figure 12: £/_ p rotates 2d partition around its diagonal. 



diagonal slice x — y — t to the 2d partition in the plane intersecting z = plane along the 
line y = — t + (1 and x = plane along the line y = —t — pz. We illustrate the action 

of U-p in Figure El We regard fl62(l as an indication that in the melting crystal picture the 
initial lattice A( p ) G (9 + has vacancies in the y-z plane at the sites: 

(0, y,z) y < N Q - pz (65) 

Initial lattice for the case of p = —2 is shown in Figure IT31 

5 D-branes in the crystal picture in the limit g s — > 

In this section we demonstrate that our proposal for D-branes in the crystal picture is 
consistent with string theory in the semiclassical limit g s — > 0. As shown in [10, 15] the 
lattice corner disappears in the limit g s — > and smooth geometry, called the limiting shape, 
emerges. The average configuration contains many molten atoms and is given by the volume 
above the 3d graph 8 

(x,y,z) = (R,U + R,V + R), U = y - x, V = z - x, R = R(U,V) (66) 

8 Thcre is also x, y, z symmetric parametrization of the limiting shape [10] 

24 




>:xx x x x xx x x x 
x x x x xxx xx xx, 

* XX XXXXX 

X X X X xxy» • • 
xxx x^» • • • • 

*-x- 



Figure 13: Initial lattice A(_ 2 ) for a D-brane with p = —2, iVo = 3. 



R(U,V) is defined as 



R(U,V) = 4 L s I 1^ log\F{U + ia,V + ip)\dadp 

4:7l 2 J JO 



(67) 



where 

F(u, v) = e~ u + e~ v + 1 (68) 

It is quite remarkable [1] that the same function F(u,v) enters the equation for the Calabi 
Yau 3-fold mirror to C 3 

Z 1 Z 2 = F{u,v) (69) 
Note that the 3d graph intersects (y, z) plane in + over the curve 



e -y + e ~ z = l, ye[0,oo] 



(70) 



Zcrystal(q) ~> 6Xp(-^J 



(71) 



In the absence of D-branes the partition function behaves as: 

Vo 

9l 

where Vo is the volume of the complement of the limiting shape ()66|) . Now let us include 
a zero- framing D-brane of the geometry (|19p. We treat this as a small perturbation to the 
melting crystal and view the additional D-brane defect in the background of the molten 
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y 



Figure 14: Area .4(iVo, 0). 

crystal. We have proposed in section 3 that in this case the melting starts from the lattice 
in + with all the sites (0, y,z), y < N g s already vacant. This implies that the deleted 
region in the presence of the D-brane Vd differs from Vq by the amount 

V D = V -g s A(N ,0) (72) 

where A(Nq, 0) stands for the area of the figure in yz-plane in + below the curve (see Figure 

M) 

e-y + e~ z = l } ye[0,N g s } (73) 

But it is exactly A(Nq, 0) that was shown in [14] to be the genus zero open string free 
energy in the A-model on C 3 in the presence of a zero-framing D-brane. This proves that 
our proposal in section 3 for the zero- framing D-branes in crystal melting picture is consistent 
with the string theory results. 

From [14] we also get evidence for the crystal interpretation of framing proposed in section 
4. As explained in [14] the genus zero open string free energy in the presence of a D-brane 
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Figure 15: Area A(N , -2). 

with framing p is given by the area A(Nq,p) of the figure in yz-plane in + below the curve 
(see Figure fTBI ) 

e -y + e- z = l, ye[0,N g s -pz] (74) 

This is consistent with our proposal in section 4 that the melting starts from the lattice 
in + with all the sites (0, y, z), y < N g s — pz already vacant. Indeed, in this case the 
deleted region would have the volume 

V D , p = V -g s A(N ,p) (75) 
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